function [output] = df(lambda)

p0 = -6;

output  = 1 - 9.8*((sech(abs(p0)/(lambda^.5)))^2)*(-abs(p0)/(2*lambda^1.5));
end
